Controlling bi-partite entanglement in multi-qubit 

systems 



Martin Plesch\ Jaroslav Novotny^, Zuzana Dzurakova^ 
Vladimir Buzek^'^ 

^Research Center for Quantum Information, Institute of Physics, Slovak Academy of 
Sciences, 845 11 Diibravska cesta 9, Bratislava, Slovakia 

^Department of Physics, FJFI CVUT, Bfehova 7,115 19 Praha l,Czech Republic 
Faculty of Informatics, Masaryk University, Botanicka 68a, 602 00 Brno, Czech 
Republic 

PACS numbers: 03.67.-a, 03.65.Bz, 89.70.+C 

Abstract. Bi-partite entanglement in multi-qubit systems cannot be shared freely. 
The rules of quantum mechanics impose bounds on how multi-qubit systems can be 
correlated. In this paper we utilize a concept of entangled graphs with weighted edges 
in order to analyze pure quantum states of multi-qubit systems. Here qubits are 
represented by vertexes of the graph while the presence of bi-partite entanglement 
is represented by an edge between corresponding vertexes. The weight of each edge 
is defined to be the entanglement between the two qubits connected by the edge, as 
measured by the concurrence. We prove that each entangled graph with entanglement 
bounded by a specific value of the concurrence can be represented by a pure multi-qubit 
state. In addition we present a logic network with 0{N'^) elementary gates that can 
be used for preparation of the weighted entangled graphs of N qubits. 
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1. Introduction 

Motivated by tlie seminal paper of Einstein, Podolsky and Rosen (EPR) pp, Schrodinger 
in his paper entitled " The present situation in quantum mechanics" j2] has introduced a 
concept of entanglement. This new type of purely-quantum mechanical correlation has 
been introduced to reflect the fact that (according to Schrodinger) ^^Maximal knowledge 
of a total system does not necessarily include total knowledge of all its parts, not even 
when these are fully separated from each other and at the moment are not influencing 
each other at all." Quantum correlations have attracted lot of attention during the 
history of quantum mechanics. Bell and Clauser et al. 4j have shown that these 
correlations violate inequalities that must be satisfied by any classical local hidden 
variable model. 

The complex phenomenon of quantum entanglement has been studied extensively 
in recent years because it represents an essential resource for quantum information 
processing (see, e.g. reference Entanglement between two qubits prepared in both 
pure and mixed states is well understood by now. In particular, necessary and sufficient 
conditions for the presence of entanglement in mixed two-qubit states have been derived 
[HI Ej and reliable measures of degree of entanglement have been introduced. Among 
others, the concurrence as introduced by Wootters et al. [SI is a very useful measure 
of entanglement since it is rather straightforward to calculate and is directly related to 
the entanglement of formation. 

Entanglement properties in multi-qubit systems are, on the other hand, still 
not completely revealed. Firstly, intrinsic multi-partite entanglement is of a totally 
different nature then a "sum" of bi-partite correlations. Secondly, unlike classical 
correlations, bi-partite entanglement cannot be shared freely among many particles [Oj. 
In particular, Coffman et al. have derived bounds on bi-partite concurrencies in three- 
qubit systems, which are referred to as CKW (Coffman-Kundu- Wootters) inequalities. 
Further investigation on entanglement sharing in multi-qubit systems have been reported 
in references [101 HH Il2l HHj- In these papers special states of multi-qubit systems 
that maximize bi-partite entanglement between selected pairs of qubits in the system 
have been presented. In addition, intrinsic multi-qubit quantum correlations have been 
analyzed (see, for instance, references [TH I15j). 

Controlling the amount of shared bi-partite entanglement in multi-qubit systems 
can be used on multi-partite communication protocols such as quantum secret sharing 
[T7j or specific multi-user teleportation schemes. 

The entanglement properties of a multi-qubit system may be represented 
mathematically in several ways. Diir for instance, has introduced entanglement 
molecules: mathematical objects representing distributions of bi-partite entanglement 
in a multi-qubit system. He has shown that given an entanglement molecule, relevant 
mixed states with the corresponding entanglement properties can be found. 

An alternative possibility for representing the entanglement relations of a multi- 
qubit system is the application of entangled graphs. The entanglement properties of a 
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system with N qubits are represented by a graph of N vertexes. The vertexes refer 
to the qubits, while the edges of the graph represent the presence of entanglement of 
the corresponding pairs of qubits. It was shown in one of our earlier papers [ini that 
for every possible graph one can find a pure state, which would be represented by that 
graph. The amount of pairwise entanglement was however not taken into account. 

In the present paper we extend the concept of entangled graphs to describe the 
amount (degree) of pairwise entanglement in the system as well. Namely, we assign a 
weight to each edge of the graph, which is equal to the amount of the entanglement 
between the corresponding pair of qubits. The entanglement is quantified in terms of a 
concurrence. 

For a given state of an N qubit system, one can obviously calculate pairwise 
entanglement, thereby constructing the appropriate graph. The inverse problem, i.e. 
finding a quantum state with entanglement properties represented by a given graph, is 
more difficult. 

In Sections 3 and 4 of the paper we will present a complete analysis of existence of 
quantum states of multi-qubit systems with entanglement properties represented by a 
given particular graph. 

For a given graph, many quantum states may be appropriate per se. The graph 
itself is not, for instance, sensitive to local operations on the qubits. On the other hand, 
there exist graphs for which no suitable state can be found. The reason behind this is 
that bi-partite entanglement cannot be shared freely: e.g. the CKW inequalities form 
an obstacle. So, for instance, we cannot have an entangled graph of three qubits such 
that each pair is maximally entangled with the value of concurrence equal to unity. In 
spite of this, a positive statement can be made. We prove in the following, that if an 
additional criterion is fulfilled, namely that the weight of each edge is bounded from 
above by a certain value, a pure state corresponding to the given graph can be found. 
This bound on the weights depends only on the number of qubits in the system. We 
also propose a constructive method, how to find these states. 

It is known that an arbitrary quantum state of N qubits can be prepared using a 
sequence of single-qubit and two-qubit operations. These operations can formally be 
represented as a quantum logic network. In general one needs to use exponentially many 
resources (counted by the number of elementary gates) to prepare a quantum state of 
qubits. 

We will show in Section 5, that for preparing a system of qubits in a state with given 
entanglement properties resulting from our consideration, less resources are needed. 
Namely, a quantum logic network composed of two- and three-qubit gates enables us 
to generate the state in argument. The number of gates building up this network is 
proportional to the number of entangled qubit pairs in the system (i.e. the edges of 
the graph). In the case of an entangled web for instance (c.f. reference II^)? when 
all vertexes of the graph are connected by edges, the number of gates necessary for a 
generation of the state is proportional to iV^. 
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2. Definitions 

2.1. Concurrence 

In this paper we will use concurrence as a measure of bi-partite entanglement. This has 
been introduced by Wootters et al. |Hj in the following way: Let us assume a two-qubit 
system prepared in a state described by the density operator p. From this operator one 
can evaluate the so-called spin-flipped operator defined as 

where ay is the Pauli matrix and a star (*) denotes the complex conjugation in the 
computational basis. Now we define the matrix 

R = pp, (2) 

and label its eigenvalues (which are all non-negative), in decreasing order Ai, A2,A3 and 
A4. The definition of the concurrence is then 

C = max|o, Y^-y^-yA^-yAlj . (3) 

This function also serves as an indicator whether the two-qubit system is separable (in 
this case C = 0), while for C > it measures the amount of bipartite entanglement 
between two qubits with a number between and 1. The larger the value of C the 
stronger the entanglement between two qubits is. 



2.2. Coffman-Kundu- Wootters inequalities 

Coffman et al. ^ have recently studied a set of three qubits, and have proved that 
the sum of the entanglement measured in terms of the squared concurrence between 
the qubits 1 and 2 and the qubits 1 and 3 is less than or equal to the entanglement 
between qubit 1 and the rest of the system, i.e. the subsystem 23. Specifically, using 
the bi-partite concurrence (jHl) the state gjk between the qubits j and k we can express 
the Coffman-Kundu- Wootters (CKW) inequality as 

C12 + C'la ^ '^1,(23) • (4) 

Coffmann, Kundu and Wootters have conjectured that a similar inequality might hold 
for an arbitrary number N of qubits prepared in a pure or mixed state. That is, one 
has 

TV 

E Cl<C]-, (5) 

where the sum on the left-hand-side is taken over all qubits except the qubit j, while 
denotes the concurrence between the qubit j and the rest of the system (denoted 

J)j2_ 

as j). The maximal value of the concurrence Cjj at the right-hand side of equation El is 
equal to unity. 
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2.3. Entangled graphs 

Let us consider a system of qubits. As already mentioned, we will represent the 
entanglement properties of the system with a weighted graph with N vertexes. Every 
qubit is identified with one of the vertexes, whereas the concurrence between a pair 
of qubits is identified with a weighted edge, connecting relevant vertexes. If a pair of 
qubits is not entangled at all, there is no edge present in the graph between the relevant 
vertexes (thus, the edge with a zero weight is equivalent to no edge). The graph itself is 
defined by the number of qubits and a set of real numbers Cjj, giving the concurrencies 
between relevant pairs of qubits. 



3. Simple examples 

The simplest example of a multi-qubit system with interesting correlation properties 
was studied in the work of Koashi et al. ^2] • These authors have studied a completely 
symmetric state of A^ qubits such that all A^(A^ — l)/2 pairs of qubits in the system 
are entangled with the same degree of entanglement. It has been shown that a state 
satisfying this condition is the so-called W^-state defined as 

\W) = \N-l), (6) 

where |A^; k) is a totally symmetric state of A^ qubits, with k qubits in the state |1) and 
all the others in the state |0). The concurrence in this case takes the value 

2 

Cmax = 'J^i (7) 

that is maximal under given conditions. 

One can easily generalize this example for other completely symmetric 
configurations (e.g. for graphs with weights equal on all edges). As proved by Koashi 
et al. [12], if the value of concurrence is larger than 21 N [see equation ((Tj)], then the 
desired state does not exist. If it is smaller than 2/N then a pure state corresponding 
to the desired entangled web reads 



1^) = Vl-a2|A^; 0) + a\N- 1) . 

The desired value of the concurrence C determines the value of a real parameter a which 
reads 



ICN C 



a 



2 y Cmax 

A more complicated two-parameter example is the case of a star-shaped entangled 
graph (see reference |IH])- In this graph a given qubit is entangled with all the other 
qubits in the system while no other qubits are entangled between themselves. In 
addition, it is assumed that the strength of the entanglement between the given qubit 
and any other qubit is the same (constant). X In the reference it has been shown, 

I This is a special case of o more general graph such that all qubits are entangled (kind of an entangled 
web 22j), but one qubit (let us denote it as the "first" qubit) is entangled with the rest of the qubits 
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that asymptotically, in the limit of large number of qubits (i.e. oo), one is able to 

find a state that saturates the CKW inequalities. Thus we are able to find a state for 
every star-shaped graph in the N ^ oo limit §. 

4. General solution 

As we have mentioned earlier, it has been conjectured, that all A^-qubit states have to 
fulfil the Coffman-Kundu-Wootters (CKW) inequalities (see equation which in the 
case when the qubit j is maximally entangled with the rest of the system reads 

Vj, T.cl,<c]j<i. (8) 

k 

Any violation of this inequality means that the corresponding entangled graph cannot 
be represented by a quantum-mechanical state. Under the assumption that all 
concurrencies Ckj in equation (jHI are mutually equal, i.e. C = Ckj, we obtain from 
the CKW inequality the bound 

~ Vn 

which is definitely not achievable. To see this we remind ourselves, that in the case of the 
entangled web (all qubits are mutually entangled) the maximal value of the concurrence 
is given by equation ((7j), which represents a bound that is much lower than the bound 
that follows from the CKW inequality. 

One may proceed either by deriving tighter CKW-type inequalities that can be 
saturated by physical states (graphs). Alternatively, one can consider only entangled 
graphs with specifically bounded weights on their edges. In what follows we will study 
this second option and will restrict the consideration to those graphs in which the 
concurrence on every edge is smaller than a certain value. We will prove that there 
exists a nonzero bound on the concurrence such that all graphs with weighted edges 
that satisfy this additional condition can be realized by pure states. 

These states are of the form 



\^) = a\A) + 1^,\B.,) (9) 



where 



1%) = (|11..0i..0,..l) + |00..1i..l,-..0)) ; (10) 
\A) ^ (I00...0) + 111...1)) . (11) 

The real positive coefficients a and satisfy a normalization condition 

2a2 + 2^7j = l. (12) 

{iJ} 

with the constant concurrence Ci, while other qubits in the system are mutually entangled as well, but 
the value of the concurrence C2 is different from Ci . 



The upper bound for bipartite entanglement given by the CKW inequalities is C < The upper 



bound for star-shaped graph is C^ax — w ~ where 6 oc ^ 
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The sums in equations Q and (fT^ are taken through all pairs i < j, i,j G (or, 
equivalently, the sums can be extended for all pairs i,j & N with the restriction 'jij = 
for j < i). The high (permutational) symmetry of the state allows us to calculate 
directly the concurrence (for details see Appendix A) 



a, = max |2 (^2a7., - E ll ~ E 7^,) , o| 



(13) 



which is valid under the condition 



« > 27„,axVAr - 2 , (14) 

where 7max = maxiji^jij). 

Let us note, that the concurrence between every pair of qubits in this rather complex 
system is expressed as an analytic function of input parameters, utilizing just a single 
condition (fT^. 

The set of non-linear equations (fT^ connect parameters of the state jij 

(the parameter a is specified by gammas via the normalization condition) with the 
concurrencies of different pairs of qubits. This set of equations is strongly coupled 
in a sense that in order to calculate one concurrence one needs to use approximately 
2A^ gammas. The task now is to invert this set of equations, i.e. to find the set 
of equation defining the gammas via the set of concurrencies that are given (these 
concurrencies do specify the character of the entangled graph). Not for every possible 
choice of concurrencies there exist parameters 'jij satisfying the normalization condition 
Ej j I'jijf < 1 and the condition (fT^ . The reason is that even though the concurrencies 
under consideration have to fulfil the CKW inequalities these inequalities are just 
necessary but sufficient condition for the existence of an entangled graph with weighted 
edges. Hence, it is also an interesting question, for which set of concurrencies one can 
find solutions of the reversed equations (fT!^ . 

We have found the solution for the parameters 7jj as functions of the concurrencies 
Cij (weights on the edges of the entangled graph) that specify the state ©, providing 
all concurrencies are smaller than a certain maximal value 

Cij < Cjnax) (15) 

where Cmax is a given constant. || . 
Theorem 1 

Every entangled graph with weighted edges that is specified by the set of concurrencies 
{Cij}, that fulfil the condition (fTHjl . can be represented by a pure state given by 
equation (jH)) 

The complete proof of this Theorem can be found in Appendix B. Here we just 
sketch how the relevant parameters 7jj can be obtained via an iteration algorithm. Let 
us start from a specific state corresponding to the situation when 

Cij Cmax 

II The upper bound for Cmax is obtained from conditions for an iteration procedure as defined in 
Appendix B. 
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for all i,j and then adjusts iteratively the parameters to fit the concurrencies. We 
can summarize the iteration process as follows: 

• After each step, all concurrencies that are evaluated for the state (jnj are greater 
than or equal to the desired set of concurrencies Cij. 

• After each step, all gammas are smaller than or equal to their values at the previous 
step; they do not change only if for a specific j the relevant concurrence is reached. 

• The iteration limit, when all gammas are zero, leads to zero concurrencies, too. 
Therefore, one has to cross the searched state during the iteration procedure (for a 
finite precision this stage can be achieved after a finite number of iteration steps) 

The existence of the state itself is proved by showing, that the iteration process 
has a proper limit. Also, to ensure the validity of the proposed process, we made a 
broad numerical test, with varying number of qubits and the strength of entanglement. 
In all tested examples that satisfied the condition (fT3j) . a very rapid convergence was 
observed, when a precision of about 10~^ of the maximal permitted concurrence was 
achieved after nine to twelve steps (changing all gammas at once). 

5. Preparation of entangled graph with weighted edges 

In the previous section we have shown that a large class of entangled graphs with 
weighted edges can be represented by a pure state 0. It is well known (see e.g. Ref. jS]) 
that any state of a multi-qubit system can be prepared with the help of a suitable logic 
network. However, in general the number of two-qubit gates in this network increases 
exponentially with the number of qubits. 

In what follows we present a quantum logic network for preparation of the state 
0, corresponding to a given weighted entangled graph. This network is very efficient 
in a sense that it uses only a quadratic number of three-partite gates with respect to the 
number of qubits (every three-qubit gate can be decomposed in at most eight two-qubit 
gates). Three ancilla qubits are needed for the procedure; these are not entangled with 
the other ones at the end of the preparation process. This keeps the fidelity of the 
preparation (in the case of error- free gates) perfect. 

5.1. Definitions 

Firstly let us introduce logic gates that will be used in our network. The first gate is 
a two-qubit operator, the well-known controlled NOT {cNOT) gate. In this gate the 
first input qubit serves as a control. The NOT operation is applied on the second qubit 
when the control qubit is in the state otherwise the second qubit does not change. 
The operator which implements this gate acts on the basis vectors of the two qubits 
under consideration as follows 



cNOT \0)^\0)^ = |0)J0)^. 
cNOT \0)^\1)^ = |0)J1)^. 



(16) 
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c7vor|i).|o). = 
c7vor|i),|i). = |i),|o)., 

where i denotes the control and j the target qubit. 

The second gate we are going to use is a three-qubit Toffoh gate T with two control 
qubits. In the case that these two control qubits are in the state |11) then the NOT 
operation is applied on the third qubit. In all other cases the Toffoli gate acts as an 
identity operator. 

The third gate we will use is also a three-qubit gate, denoted as R{a). Here one 
qubit will serve as a control. When this control qubit is in the state |1) then a specific 
"rotation" in the two-dimensional subspace of the Hilbert space of the two target qubits 
will be apphed. This rotation acts on the two target qubits as follows: 

R (a) |00) = (l - {aff^^ |00) - a |11) ; (17) 

R{a)\n) = a |00) + (l - (af) 1 11) ; 
R{a) |01) = |01) ; 
R{a) |10) = |10) , 

where ck is a parameter of the rotation. 

The last gate that will be used in our network is a three-qubit gate, which will be 
denoted by A. Again, one qubit serves as a control. The following transformation is 
applied to the two remaining qubits when the control qubit is in the state |1): 



A |00) = k+ |00) - k+ |01) - A;_ |10) + A;_ |11) ; (18) 

A |01) = A;_ |00) + k+ |01) - A;_ |10) + k+ |11) ; 

A |10) = k+ |00) + k- |01) + k+ |10) - k^ |11) ; 

A 111) = - fc_ |00) - fc_ |01) + k+ |10) + k+ 111) , 

where we have used a short-hand notation k± = \\J^ ± This operation will be used 
in our network even number of times, so only the effects of the operation A^ will appear 
at the end. The operation A^ acts in a simpler and understandable way 

A^ |00) = - |01) ; (19) 

A'm = 111) ; 

|10) = |00) ; 

^2 111) = |10) . 

5.2. Initial state of the N qubits 

In order to prepare an entangled graph with vertexes, i.e. a specific A^-qubit state, we 
will need three additional ancilla qubits. The ancilla is initially prepared in the product 
state |1)|0)|0) and is completely factorized from the other, graph qubits. These graph 
state are initially prepared in the generalized Greenberger-Horne-Zeilinger (GHZ) state 
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Figure 1. A schematic description of the logic network corresponding to the first stage 
of the preparation procedure. The rotation R{a) is apphed on every pair of the target 
qubits (all original qubits, except the ancilla qubits), with a defined by the parameters 
7.y that specify the state corresponding to a given entangled graph. During this first 
stage of the preparation procedure approximately iV^ elementary gates are used. 

= — |0)®^)/v^. Thus the input state of the quantum logic network under 

consideration reads 

1^)/ |0)^+2 10)^+3 = ^ (111-1) - I00...0)) |0);v+2 |0)iv+3 -(20) 

In what follows we will specify gates in the network with three indices, where the 
first index specifies the control qubit (or the first two qubits in the case of the Toffoli 
gate) and remaining index(es) determine(s) the target qubit (s) of the operation. In 
addition, if there will be some action or control applied on the ancillas, we will denote 
their relevant indexes as iV + i, where z = 1, 2, 3 is the position of the ancilla qubit. 

5.3. The network 

The action of the network can be divided into two main stages. In the first stage an 
entangled state of the graph qubits and the ancilla is created. This state contains state 
vectors that are essentially the same as those in the desired state ©• In the second 
stage of the preparation procedure the ancilla becomes factorized from the graph, which 
in turn is prepared in the state 

During the first stage of the preparation procedure we will apply the rotation R {aij) 
to each pair i ^ j from target qubits with the control on the first ancilla qubit (see 
Figure^). After each R (ajj)-gate the Toffoli gate with the control on i and j qubits acts 
on the first ancilla qubit. This procedure is repeated (^^^-times, for all indices i ^ j. 
During each rotation, a fraction (that is specified by the amplitude aij) of the state 

% To generate a GHZ state one can start with a product state of N qubits with the first qubit in the 
state (|1) — |0))/V2 while all other qubits are in the statejO). Then one applies a cNOT gate to every 
qubit except the first one with the control on the first qubit. So, one needs only — 1 two-qubit gates 
to prepare the input state. 
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vector - I00...0)) is transformed into the state (|11..0i..0j..l) + |00..1i..lj..O)), 

whereas the aheady transformed part of the state (|11..0fc..0;..l) + |00..1fc..li..O)) is left 
unchanged. 

Thus, after a few steps the state given by equation (j^Uj) is transformed into 

m = (|11..0...0,..l) + |00..1...1,..0)) |0) |0) |0) (21) 

{id} 

+ Z}(|11...1)-|00...0))|1)|0)|0), 

where the tilde indicates that the sum is taken over all pairs of qubits that have already 
been involved in the transformation. The corresponding amplitudes Cij are given by 
the relation 

E., = ^il-2Y.ElX (22) 

\ {k,l} / 

and D is given by the normalization condition. 

We have also to specify the parameters aij for each rotation. These parameters 
are related to amplitudes 7,^ that specify the desired state of the entangled graph given 
by equation 0, which we want to generate. Comparing the states ^ and we 
see, that for a successful generation of the state of the graph we need Eij = •jij. Using 
equation we can write 



Oiij 



V2 1-2^:7.'^ l^r (23) 

V J 

After performing transformations on all pairs of target qubits the resulting state 
has the form 

\^) = Y: 7^,\B^,) |0) |0) |0) + a \a) |1) |0) |0) . (24) 

where the state vector \Bij) is given by equation (fTUj) and 

\A)^ (|11...1)-|00...0)) . (25) 

We see that the component states \Bij) and \A) in equation are essentially the 
same as those of the desired entangled graph (see equation Q). Now we will use the 
first ancilla qubit the last time before disentangling it from the rest of the system. We 
will apply the specific controlled rotation on an arbitrary qubit of the graph with the 
control being the first qubit of the ancilla. The rotation itself is described by the operator 
—(7z (a Pauli matrix). This controlled rotation applied on the state performs the 
transformation \A) — > \A), while the state \Bij) remains unchanged. 

We see that at this stage the two desired components a\A) and J2lij\Bij) of the 
graph state ^ are generated, but they are entangled with the first ancilla qubit. The 
second stage of the preparation procedure is designed so that the ancilla is disentangled 
from the graph, while the graph is left in the state ©. To disentangle the first ancilla 
qubit from the rest of the system we will use the other two ancilla qubits. In order to 



Controlling bi-partite entanglement in multi-qubit systems 12 



— • 


• Q 






B 






































IP 


© 














w 






















1 


i 


V 


1 









Figure 2. This part of the network helps to disentangle the ancilla qubits from 
the original qubits. In this case, the rotation A is applied only in the case, when 
two neighboring qubits are not equal. Only N elementary gates, in the order of the 
magnitude, are used. 

perform this disentanglement we have to find a network that will discriminate between 
two graph states \A) and 

Let us analyze in more detail the state J2lij\Bij). In this state two or four 
neighboring qubits are in different states. This is in contrary to the state \A), in which 
all qubits are in the same state. The discrimination of the two states can so be performed 
by applying the cNOT gate acting always on two neighbouring qubits. If the target is 
in the state |0) then after the action of the cNOT gate the two qubits that are involved 
in the action of the gate are in the same state. On the other hand, if it is the two 
qubits do differ. From here it follows that we can use the target qubit as a control for 
another gate, which changes its targets if and only if the two graph qubits do differ. 

Let us utilize for this purpose the gate A which will act on the last two ancillas (see 
Figure El). We apply the cNOT gate on two qubits from graph qubits and then we 
apply the A-gate controlled by the target of the cNOT, acting on the last two ancillas. 
After that, we apply again the cNOT gate on the same two qubits as before: This 
operation will bring all qubits into the original state (since cNOT^ = I) and the only 
effect of this particular procedure is a rotation of the state of the last two ancillas. This 
rotation will take place only in the case when the two "tested" qubits were in a different 
state. 

Than we repeat the same procedure for each pair of the first neighboring qubits 
of the graph. After this, the A gate acted either 2 or 4 times on ancilla qubits, that are 
entangled with the state J^lijlBij) of the graph. On the other hand those ancillas that 
are entangled with the state 1^4) are not changed. 

The reason for using A-gate, the "square root" of the operation (|T!l|l . now becomes 
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Figure 3. A schematic description of the entire logic network for the preparation of 
the entangled graph in the state First we use the part with the rotation R{a), then 
we correct a sign with the help of a cNOT gate. Finally, we disentangle the ancilla 
qubits with the help of the A and A^^ rotations. The desired state is prepared with 
the help of (of the order of magnitude) iV^ elementary gates. 

clear: the gate is acting once or twice and the state |0) |0) of the last two ancillas is 
changed either to — |0) |1) or to |1) On the other hand in the case when all target 
qubits are equal, A will not act at all and the resulting state of the last two ancillas will 
be unchanged, thus |0) |0). If the cNOT gate between the last (control) and the first 
(target) ancilla qubit is applied, then the first ancilla will be changed to the state |0) 
and it becomes disentangled from the rest of the system. Now all the work is almost 
done, the only thing we have to do is to disentangle the two remaining ancilla qubits. 
For this we will simply run the procedure for all neighboring pairs of qubits as described 
above, but with the gate instead of the gate A. This will change the state of the 
last two ancilla qubits back to the original state |0) |0) and will finally disentangle the 
ancilla from the system. That means that the desired state l^') of the graph qubits 
is disentangled from the ancilla and the entangled graph is prepared in the state (jnj. 

Finally, let us summarize the preparation procedure for the entangled graph given 
by the state (jH))- As shown in Figure El first we apply the rotations R on all pairs of 
qubits of the graph, i.e. 

Gi = (ccT^)7v+i,iv n {Ti,j,N+iR (aij) N+i,i,j) ' (26) 

where the subscripts for each operation define the position of qubits, where operation 
takes its action. Angles of rotations a^j are defined by equation and ca^ stands for 
the controlled sigma gate applied on the first ancilla as a control and one of the graph 
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qubits as a target. At this stage we will use roughly A^^ bipartite gates. The second 
stage of the preparation corresponds to disentangling the first ancilla qubit from the 
graph qubits 

N 

G2 = cN0T^+3,N+i n {cNOT,,,+iA+,,^+2,N+3cNOTi,i+i) . (27) 

The last stage of the preparation process is responsible for disentanglement of the last 
two ancilla qubits from the graph qubits, i.e. 

N 

G3= U (cArOT,,,+iA-^+i,;v+2,Jv+3cArOT,,,+i) . (28) 

In the last two equations the indices z + 1 for the gates are taken implicitly as modulo 
A^. Finally we can represent the action of the whole logic network as 

|^)i.|l),V+l|0)7V+2|0)^+3 = G'3.G'2.G'i.|^)/|l)^+i|0),V+2|0)^+3, (29) 

where is the desired state © of the entangled graph with weighted edges. 
6. Conclusion 

In this paper, we have introduced a concept of entangled graphs with weighted edges. 
Using simple examples we have shown, that sharing of bipartite entanglement is a 
complicated phenomenon and that the Coffman-Kundu-Wootters inequalities [U] are 
only a necessary condition for an existence of states with given entanglement properties. 

We have proved that a whole class of entangled graphs, where the concurrence 
between an arbitrary pair of qubits (vertexes) is weaker than a certain value can be 
realized by a state of A^ qubits. Moreover, we have proposed a logic network for 
preparation of the states corresponding to this entangled graphs. The network is 
composed of a number of elementary quantum gates that grows quadratically with the 
number of vertexes (qubits) in the graph. 
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Appendix A. Concurrence in entangled graphs 

In what follows we will evaluate the concurrence between an arbitrary pair of qubits of 
a system in the state Q, i.e. 

|vl/) = a (I00...00) + |ll...ll)) + 5: 7., (|11).,|00...00)^+ |00).,|ll...ll)^) , (A.l) 

{iJ} 

where real positive amplitudes a and satisfy the normalization condition 

+ 2 E 75 = 1. (A.2) 
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The sum in equations ()A.1|) and ()A.2|) is taken through all pairs i ^ j, i,j G N, so 
{i,j} = {j, i} and thus •jij = for z < j. The special form of the state ()A.1|) leads to a 
rather compact density matrix for an arbitrary two-qubit operator that is obtained by 
tracing over the rest of the graph qubits: 

/AO F \ 
B E 



E B 
\ F A / 

where we have used the notation 

,2 



(A.3) 



A 



{fc,/} 



(A.4) 



2 E 7fci7jfc ; 

k 



E 
F 

All sums in equations ()A.4|) are running through free parameter (s) k (and I), whereas 
i and j denote a specific pair of qubits in the graph. In addition, the condition 
i k I j has to be fulfilled. 

The convenient form of the matrix ()A.3|) allows us to calculate square roots of the 
eigenvalues of the matrix R given by equation (j2I): 



Ai 

A2 

As 
A4 



A + F 
A-F 
B + E 
B-E 



(A.5) 



Because the coefficients A, B, E, F are positive, the only candidates for the largest 
eigenvalue are Ai and A3. Let us further define 



7max = max^'jij] 



Using the condition 

«>27n 



(A.6) 



(A.7) 



we find Ai > A3 and the general expression for the concurrence associated with edges of 
the entangled graph prepared in the state (jA.ljl reads 



Cij = max ^ 2 ( 2a7j 



E 

{k,i} 



E 7fci 







(A.8) 
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Appendix B. Proof of Theorem 1: Iterative procedure 

In order to prove Theorem 1 we first label the set of concurrencies that determine a 
given entangled graph by Cjj. We will use a bold C in order to distinguish these 
concurrencies from any intermediate concurrencies, obtained by searching for the state 
of the entangled graph. 

We will start the iteration procedure with an initial state of the entangled graph 
given by equation 0. The amplitudes are specified by the relation 

, ^ (B.l) 

that is, the initial state is completely permutationally symmetric. The parameter A is 
defined as 

^A{N - 2Y + 2N{N - I) - 2{N - 2)) 

^ = N{N^) • ^^-2^ 

The corresponding bi-partite concurrencies can be evaluated straightforwardly and they 
read: 



C§^ = C^^ = 2{a^hS^-2iN-2)(stS^r) 



VQN^ - 18N +16-2N + 4: 
= N{N - 1) • ^^-^^ 

We remind ourselves that the parameters a^^^ and 71°^ are mutually related via the 
normalization condition ()12|) . therefore a is always implicitly defined by 7jj. It is also 
clear that for the state under consideration the condition (fT^ is fulfilled as well. 

Before we describe the iteration procedure itself we introduce the following notation: 
we enumerate all pairs of qubits in the entangled graph. All pairs of qubits (i.e. the 
edges of the graph) are listed in the set of pairs just once. At each iteration step one 
parameter 7^^ for a selected pair of indices {k, 1} is changed, whereas all others gammas 
will stay unchanged. Let us now suppose, that the n-th step of the iteration is done and 
both conditions (fTH|l and (fT^ are still fulfilled. Moreover a^"'\'yl^'^ are positive. Hence 
we find 

> ; (B.4) 

a(")>2v/iV^7il; (B.5) 
< < 1 < 7^;^ < 1 , (B.6) 

for all pairs of indices i,j. The parameter 7^2x is defined in the same way as in 
equation ()A.6|) . i.e. 

7tl = max(7l;^ (B.7) 

In the next iteration step we take a pair of qubits (i.e., the edge), that follows after 
the pair, which was selected in the previous iteration step n. Let us denote this pair with 
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indices {«, j}- Then, in the {n + 1) — st iteration step, we will change the parameters of 
the state in a following way: 

2 



Tir = . ; (B.8) 



Tjiri) _|_ T/(n) 



where 



1/2 



, . 1/2 



(B.IO) 



(B.ll) 



All other gammas remain unchanged at this iteration step. Following the conditions 
()B.4|) and ()B.5|1 this iteration step is well defined. Now we will discuss several important 
properties of the iteration process: 

(1) a*^""*"^^ and 7^"'*'"^'' are solutions of the equation 

^^"-^^^'^ = + 7 (C.. - (B.12) 



and thus according to equation (fT!?|) 



C^^^ = max 2 2a("+Mr^) - E (Tif'V " E (7^'^)' | ' 

I \ {k,i} {k,j} 

= max{Qj,0} = Q^. (B.13) 

(2) a^""*"!) and 'JiJ'^^^ fulfil the normalization condition ()12|) . 

(3) 7^"^"'^'* and a^'^"^"^^ are positive and satisfy the relations 

<7;r^<7j?; (B.14) 
< < 1. (B.15) 

(4) From equations ()B.14j) and ()B.15|1 it follows that 

> ^(n) > 2v/iV^^W^ > 2v/iV^7tt'^- (B.16) 
Therefore the condition ()B.5jl is valid also for the (n + 1) — iteration step. 

(5) Let us now show, how will particular concurrencies change in this single iteration 
step. For k, I ^ i,j we find 

^ = 2 f2a(-Sir) - E (71^^^- E (7r^Tl ^^-1^) 

\ {k,m} {l,m} / 

=2(2a(-si?-E(7S)'-E(7;^y 

>2f2a(")7i?- E (7^2)'- E (7tT 

\ {k,m} 

_ 
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and for k = i 

{j,m} {^™} 

> 2 (2aH7i? - E (7tT - E (7.t^ 

(n) 



= 0; 

The same is valid also for = j 



Thus we have shown, that after this iteration step the concurrence for fixed j (i.e. 

-y(n+l) _ f-^ 



for the given edge) will be C^"''^^^ = C^- and all other concurrencies of the entangled 



graph will become larger. Thus, the condition for all i,j C}"^^'' > Cij will be fulfilled. 
Therefore, the state defined by equation (jH)) with the parameters 7^"'*'^'' can be used for 
the next {n + 2)-nd iteration step. 

Therefore, the whole iteration is well defined and we will obtain an infinite sequence 
of parameters (a'-"'''} and |7,-!^^} for each pair of indices ij (i.e. for each edge of 

L J n=0 L J n=0 

the entangled graph). All sequences are monotonous and are bounded, and therefore 
they have proper limits. Let us denote these limits as a and jij 

a = lim a^") ae (0,1) (B.19) 

n^oo 

7^.= 1H^7;;^ ^ 7.,G(0,1). (B.20) 



n— >oo 



Now we will choose and fix one pair of indices i,j and we will show, that 

lim C\f = Q,. (B.21) 

First we define a sequence {k{n)}'^^Q in a following way: fc(l) = p, where p is a rank of 
{i,j} in the order of pairs of indices, and k(n) = p + "^(^"^) , Then 

= C,,. (B.22) 
The equation ()B.21|) is equivalent to the definition 

(V£ G R, e > 0) (3no e N) (Vn G N, n > no) (l^.^"^ - Cij\ < e) . (B.23) 

Let us choose and fix the small parameter e. Our task is to find no, that will have the 
property ()B.23|1 . Because all sequences ja*^"")! ^ and |7i"''| ^ have a proper limit, 
they are Cauchy sequences and therefore 



(VrGR,r>0)(3moGN)(Vn,mGN,n,m>mo)(V{fc,/}) I („) ) , (B.24) 



where 



m{N - 1) ■ 



(B.25) 



Controlling bi-partite entanglement in multi-quhit systems 



19 



For this r there exists such mo, that the property ()B.24|) is fulfilled and we can define 

uq as 

no = k{mo) > mo (B.26) 

(n+i) ^(n)| for n + 1 > no and 

^ Q. J. XiV/ ACLOU ^WAJ-V^IUAWAJ. A J-AV^CLXJ-O , UAJ.CLU UAJ.\^ yi L | Xy itCratiOU StCp 

didn't change 

difference under consideration, either 



Further we will calculate the difference |C/" — C„ 
n + 1 ^ {^('^)}^o- '^^^ ^^^^ condition means, that the (n + 1 

From equations ()B.17j) and ()B.18j) we obtain two options for the 



^(n+l) ^{n) 



4 



a 



(n+l) 



a 



in) 



(n) 
lij 



< 4r < 8r . 



(B.27) 



or 



I a 



(n+l) 



(n)| 



+ 2 



(n) 



(n+l) (n) 
Til - Til 



(n+l) . (n) 
7ii + 7ii 



< 4 

< 8r, 

where 7^^;"''' is the parameter, which was changed in the (n + 1) — st iteration step 



Finally, we can say for n > no, if n G {^(^)}^o' then |Ci"^ 



C,J = 0. In the 



opposite case there exists such u G Nq, that 
n G (A;(mo + u), k{mo + u + 1)). 

Thus 



Cij 



< 



\N{N-1)t 



.(n-1) 



a 



(n-2)^ 



+ 



(B.28) 



/^(/c(mo+M)) ^, 



But then it must stand 



lim a 



(n) 



Jim2|2aW7i;)-E(7i?)^-E(7i?^^ 

{k,j} 



2|2«7.,-E7.'.-E7^, • 

All other conditions remain fulfilled in the limit form as well. Because this property is 
valid for all pairs of indices, we have found the parameters jij, that define the state Q 
which corresponds to a given entangled graph. 
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